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We study the existence of D-brane bound states at threshold in Type II string theories. 
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^ ■ In a number of situations, we can reduce the question of existence to quadrature, and the 
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study of a particular limit of the propagator for the system of D-branes. This involves 
a derivation of an index theorem for a family of non-Fredholm operators. In support of 
the conjectured relation between compactified eleven-dimensional supergravity and Type 
IIA string theory, we show that a bound state exists for two coincident zero-branes. This 
result also provides support for the conjectured description of M-theory as a matrix model. 
In addition, we provide further evidence that there are no BPS bound states for two and 
three-branes twice wrapped on Calabi-Yau vanishing cycles. 
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1. Introduction 

Remarkable progress by Polchinski in describing the solitons of Type II string theory 
has provided the means by which many conjectured dualities involving string theories and 
M-theory can be stringently tested [Q. The low-energy dynamics of coincident D-branes 
has been described by Witten 0, who reduced the question of finding BPS bound states 
to one of studying the vacuum structure of various supersymmetric Yang-Mills theories. 
In simple cases, the BPS mass formula forbids the decay of a charged particle saturating 
the mass bound; hence, ensuring stability. However, there are a number of situations in 
which a particle is required that is only marginally stable against decay. Showing the 
existence of such particles, with energies at the decay threshold, is the goal of this paper. 
A similiar problem arose for finite SU(2) N=2 Yang-Mills theory, studied in and ||, 
where certain dyon bound states at threshold were shown to exist. The situations we shall 
presently study are significantly more difficult because the Hamiltonians are not as well 
behaved, and gauge invariance provides an added complexity. 

Let us briefly recall the low-energy dynamics of coincident Dirichlet p-branes, de- 
scribed in 0. The world volume theory is the dimensional reduction of ten-dimensional 
N=l Yang-Mills to the p + 1-dimensional world volume of the brane. For a single brane, 
the gauge theory is abelian, and the dynamics therefore trivial in the infrared. For N 
coincident branes, the gauge symmetry is enhanced to U(N) rather than U(l) . After 
factoring out a U(l) corresponding to the center of mass motion, the existence of a bound 
state requires that the remaining SU(N) p + 1-dimensional Yang-Mills theory possess a 
normalizable supersymmetric vacuum. The bosonic potential for this model generally has 
flat directions, and so we encounter the problem of bound states at threshold. If a bound 
state is required by a conjectured duality, there is a consistency check, described by Sen 
|J, that can sometimes be performed. In favorable cases, one might be able to further 
compactify one direction of the superstring theory. If a bound state exists prior to com- 
pactification, it should give rise to BPS states in the further compactified theory which, 
for appropriate choices of momentum along the circle, are no longer marginally bound. 
The existence of these states can then be analyzed with more conventional techniques. Of 
course, for this consistency check, there have to be enough remaining uncompactified di- 
rections so that problems with infra-red divergences do not arise. More generally, however, 
the question of bound states at threshold must be addressed. Note that a normalizable 
state for a theory in a compact space generally does not remain normalizable when the 
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volume is taken to infinity. The spectrum can and often does change discontinously, and 
showing the existence of the bound state in the non-compact situation requires a separate 
analysis. 

In a similar spirit, we can arrive at descriptions of the effective dynamics of p-branes 
multiply-wrapped on supersymmetric cycles of a compactification space. In the case of p- 
branes wrapped on p-cycles, the resulting description of the low-energy dynamics is some 
flavor of quantum mechanics, although not generally just a supersymmetric gauge theory. 
Our aim in this paper is to address the fundamental issue - the existence of flat directions 
in the potential - which arises in studying binding in these situations. This analysis 
generalizes the discussion in ||, where we argued for the existence of a marginal bound 
state of a zero-brane and a four-brane, to the case where the gauge group is non-abelian. 
We shall see that there are very subtle issues that arise as a result of this complication. 

The most exciting reason for studying this question is, however, the remarkable con- 
jecture that M-theory may be described in terms of zero-brane dynamics in the limit where 
the number of branes goes to infinity 0. This conjecture is, in part, founded on previous 
work studying the relation between supermembranes, and the iV — > oo limit of type IIA 
zero-brane quantum mechanics HQ. In order for the M(atrix) model to have a chance at 
describing M theory, we need to be able to find states in the quantum mechanics which 
correspond to the gravitons of eleven-dimensional supergravity. The bound state that we 
shall find is precisely one of these particles. In the process of showing that such a bound 
state exists, we will provide a detailed study of the behavior of the propagator for the two 
zero-brane system when the zero- branes are far apart. There are a number of complications 
that make this analysis quite subtle. During the lengthy course of our investigation (which 
pre-dates the M(atrix) model), a number of germane papers have appeared. Among these 
papers have been interesting discussions of zero-brane scattering in various approximations 
and more recently, an exciting extension of the original matrix model conjecture 



to the case of finite N ||12|| . There has also been an explicit argument showing that there 



are no normalizable ground states in a particular simplified matrix model |13[], a heuristic 
attempt to argue for the existence of zero-brane bound states |jH[ , and a recent paper 
which has some overlap with our results [ITH . 



In the following section, we consider the case of p = 0. We describe a seven parameter 
family of theories, which are the primary focus of this discussion. This family of theories is 
derived fundamentally from the quantum mechanics describing the zero-brane in type IIA 
string theory by adding mass deformations. These parameters allow us to 'flow down' from 
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ten dimensions to models that correspond to the reduction of N=l Yang-Mills in lower 
dimensions by taking various mass terms to infinity. We discuss general features of these 
models, including the various physical scenarios in which they arise. Our approach to the 
question of counting bound states is described in section three. There, we argue that the L 2 
index for this class of supersymmetric quantum mechanical Yang-Mills theories is actually 
computable. This involves a discussion of L 2 index theory for non-Fredholm operators, 
which is an area of mathematics that is relatively unexplored. In section four, we study the 
question of two-particle binding in these models, and we derive a formula for the principal 
contribution to the index. The final section is a study of the two-particle propagator in the 
limit where the two particles are far apart. With this analysis, we can compute a subtle 
additional contribution to the index. The way in which this contribution arises involves 
some rather surprising cancellations. In the class of models that we investigate, we find that 
only the case which corresponds to the reduction of supersymmetric Yang-Mills from ten 
dimensions can have a unique bound state. This answers, in large part, the question of why 
the large iV limit of the reduced ten-dimensional Yang-Mills theory should be distinguished 
from the large iV limit of reductions of lower-dimensional Yang-Mills theories. 



2. Quantum Mechanical Gauge Theory 

2.1. General Comments 

Let us begin by considering models that arise from reducing supersymmetric d+1- 



dimensional SU(N) Yang-Mills to quantum mechanics; see, for instance, [16[ for the first 
discussion of quantum mechanical gauge theories, or perhaps ||17|| . Whether the Yang-Mills 
theory contains additional matter multiplets does not significantly change the following 
discussion; so, for simplicity, we shall assume no additional matter. On reducing the 
connection A M , we obtain scalar coordinates x % where i = 1, d which take values in the 
adjoint representation of the gauge group. We introduce canonical momenta obeying, 

where the subscript A is a gauge index. With the generators T A for the adjoint represen- 
tation normalized so that Tr(T A T B ) = N5 AB , the Hamiltonian for the system takes the 
general form, 

H = 2^Tr(pV) + V(x) + H F . (2.1) 



The bosonic potential V(x) is polynomial in x, and generally has flat directions. The term 
Hp is quadratic in the fermions and linear in x. Specific examples will be studied in the 
following subsection. The Aq equation of motion gives a set of constraints, Ca, which must 
vanish on physical states by Gauss' law. The constraints obey the algebra, 

[C A ,C B }=ifABcC c , (2.2) 

where Jabc are the structure constants. The constraints further obey the commutation 
relations [Ca,H] = [Ca,Q] = 0, where Q is a supersymmetry generator. The supersym- 
metry algebra closes on the Hamiltonian if the constraints are set to zero. An iV-particle 
BPS bound state corresponds to a normalizable, gauge-invariant ground state for this 
supersymmetric system. 

Without detailed computation, what might we infer about the structure of the ground 
state? Away from the flat points, the wave function for the ground state will decay 
exponentially. The only interesting asymptotic behavior is expected near points where the 
potential is small. A preliminary comment about the structure of the flat directions is 
in order: for gauge group SU(N), there are d c = (d — 1)(N — 1) + (A 2 — 1) commuting 
directions around a flat point, and d a = (d — 1)(A 2 — N) non-commuting directions. Let 
us consider the structure of the potential in the neighborhood of a flat point. As we 
shall subsequently describe in detail, the potential can be approximated by V ~ — ^r 2 |u| 2 , 
where v parametrizes the transverse directions, and r is a radial coordinate for the flat 
directions. The Hamiltonian is then essentially a set of bosonic and fermionic harmonic 
oscillators for the transverse directions, and a free Laplacian along the flat directions. 
The frequency for oscillation along the massive directions depends on r. This observation 
provides one way of seeing that there are no scattering states in the spectrum of the 



bosonic Hamiltonian for these models, as discussed in \TB\. Somewhat surprisingly, the 
spectrum of the bosonic models only contains discrete states. To construct a scattering 
state along the flat direction, one would want to put the transverse harmonic oscillators 
into their ground states; however, the zero point energy of the oscillators increases with r, 
essentially forbidding finite energy scattering states. The same argument does not apply 
to the supersymmetric case, since the ground state energy for the additional fermions now 
cancels the zero point energy from the bosons, as required by supersymmetry. If this were 
not the case, the subtleties in counting zero-brane bound states would not exist! 

In a first approximation for large r, any zero-energy wavefunction, ij)(x), roughly takes 
a product form corresponding to placing the transverse oscillators into the ground state, 
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ijj(x) ~ g(r, 0)e~ r \ v \ / 2 , where 9 are angular variables for the flat directions. The leading 
dependence of g(r, 9) on r is believed to be power law decay for large r. Acting with the 
Hamiltonian for the massive directions on this wavefunction yields zero, since the zero 
point energies of the bosons and fermions cancel. We can now explain the key difficulty in 
studying the approximate asymptotic wavefunction: can the decay exponent be accurately 
estimated? 

We note that this issue is critical, and cannot be resolved by simple approximations 
of the asymptotic behavior. For instance, even in this approximation, the function g(r, 9) 
is not simply the solution of a free Laplacian for the <i c -dimensional space of flat directions 
since the Laplacian, which for the radial coordinate is given by, 

*r = -~r^, (2.3) 

also acts on the harmonic oscillator component of the wavefunction. Actually, it is unlikely 
that the decay exponent can be accurately estimated without at least including the first 
excited mode for the massive direction into the approximation. We should also note that 
showing that the decay is fast enough to ensure normalizability is only a first step toward 
showing that a bound state exists. The structure of the wavefunction would need to be 
studied at small r where the non-abelian degrees of freedom are important. Currently, the 
only practical approach is to develop an appropriate index theory for the problem. As a 
final comment, note that the power law behavior of the asymptotic ground state wave- 
function is a consequence of the lack of a mass gap in the spectrum. The supersymmetric 
theory contains a continuum of states which descend to zero energy, thanks to the existence 
of the fiat direction || . 

2.2. A Family of Models 

We now turn to the models of primary interest to us. Let us recall that strongly 
coupled Type IIA string theory in ten dimensions has a conjectured dual description as 
weakly coupled eleven-dimensional supergravity compactified on an 5" 1 [ffSLOT. To match 
the Kaluza-Klein spectrum of the compactified supergravity theory, Type IIA string theory 
requires electrically charged particle states. The Dirichlet zero-branes, which carry RR 
charge, seem to be the only candidates. Since there is a single Kaluza-Klein mode for each 
choice of momentum along the circle direction, we desire a single D-brane bound state for 
each N . Proving this conjecture was our original motivation for studying these theories. 
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Actually, Sen has argued in that if a unique bound state exists in the quantum 
mechanics describing N zero-branes, then the spectrum of ultra-short multiplets in the 
toroidally compactified type II string agrees with the spectrum predicted by U-duality. 
The world- volume theory for the D-particle is given by the dimensional reduction of N=l 
9 + 1-dimensional Yang-Mills to quantum mechanics. A Majorana-Weyl spinor in 9 + 1 
dimensions has 16 real components, which means that the resulting quantum mechanical 
theory has N=16 supersymmetry. Let 7^,0 be a real representation of the SO (9) Clifford 
algebra with i = 1, 9 and a = 1, 16. These Clifford matrices satisfy 

{ 7 \y} = 2S ij . 

After reduction, the Hamiltonian for this system takes the form, 

H = 2^(PV) - ^ ^IV([x 1 ,^] 2 ) - ^T>(V>rV, V]), ( 2 - 4 ) 

ij 

where the real fermions ipAa obey: 

{tpAa^Bp} = SABS a p- (2.5) 

The Hilbert space is then composed of spinors on which the quantized fermions act as 
elements of a Clifford algebra. The spinor wavefunctions contain an extremely large number 
of components, even for small N, which makes an explicit construction of the zero energy 
bound state wavefunction at best difficult.S The supersymmetry algebra takes the form, 

{Qa, Q(3} = 26**11 + 2^ ap x\C A , (2.6) 

where, 

while the constraint, 

C = -i[x\p 1 } - -[ip a ,^a], 

or explicitly, 

C A = fABc(x B Pc ~ -^Ba^Ca)- (2-7) 

1 However, the existence of the nice Spin(9) flavor symmetry might bring an explicit construc- 
tion of the ground state wavefunction within the realm of possibility. We leave the attempt to 
construct the explicit solution to braver souls. 
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The constraint takes exactly the form assumed in the previous discussion. It is natural 
to call this a nine-dimensional model, although it is quantum mechanics, since there are 
nine bosonic variables in the adjoint of SU(N), and the model is the reduction of a ten- 
dimensional theory. The flavor symmetry is clearly Spin(9). Note that there is a nice 
correlation between fermion number and flavor representation that is worth mentioning 
at this point. The correlation essentially follows from spin-statistics in ten dimensions: 
fermionic states in the Hilbert space transform under spinor representations of the flavor 
group, while bosonic states appear in representations of SO (9). If the ground state is 
unique, it must therefore be bosonic. There are similar relations for the other models that 
we shall soon discuss. 

Let us consider what sort of deformations are possible in this theory. We would like 
to add mass terms to compactify some of the bosonic variables and effectively reduce the 
dimension, but we will also require that the supersymmetry algebra maintain its nice struc- 
ture. In particular, we shall not consider deformations which introduce additional terms 
into the right hand side of the supersymmetry algebra (2.6), which are linear in momenta. 
The mass deformations that we shall describe correspond, in special cases, to breaking N=4 
Yang-Mills in four-dimensions to N=2 or N=l by giving masses to various chiral fields in 
the adjoint representation, and reducing the corresponding model to quantum mechanics. 

To describe the allowed deformations, choose a real supersymmetry generator, Q = 
Q a . The generator can be split into terms involving momenta, and terms independent of 
momenta. Those depending on momenta can be expressed, schematically, as X\p\, where 
A* is a real fermion, and i runs from 1 to 9. This leaves us with seven real fermions, uj 3 a in 
the adjoint of the gauge group, unpaired with a momentum operator, but each appearing 
in Q paired with an operator, f A , quadratic in the coordinates. The supersymmetry 
generator is then roughly, 

q~\ aP \ ■ J A f A • .... 

The seven fermions, uj- ? , then represent our deformation degrees of freedom. We can add 
any reasonable operator to f J , independent of the momenta, and not generate a new 
term linear in momenta in the expression for {Q, Q}. There are many interesting possible 
deformations that preserve at least one supersymmetry. Some deformations can give quite 
exotic classical minima of the resulting bosonic potential. This is a topic that merits 
further investigation. As a special prosaic case, we could add the perturbation mx 1 to one 
of the f J \ which would lift some of the flat directions. This is the family of deformations 
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to which we shall restrict our discussion. More explicitly, consider a term / J which squares 
to give the term in the potential, |/ J | 2 . Adding the term mx l to f j changes the potential 
to I/- 7 + mx l \ 2 . Taking m — > oo then effectively decouples x 1 from the model. In this way, 
we generate a seven parameter family of models which depend on the values of the allowed 
masses for seven of the coordinates. Note that taking all masses to infinity leaves us with 
a two-dimensional model, and further compactification is not possible without introducing 
additional terms linear in the momenta into the supersymmetry algebra. 

There are two cases of particular interest: the three and five-dimensional models. 
These models correspond to the reduction of N=l Yang-Mills from four and six dimensions, 
respectively. For completeness and to fix annoying normalizations, we shall describe the 
Hamiltonian and supersymmetry algebra for both models explicitly. 

In the three-dimensional case, the Hamiltonian is given by, 

H = 2^' &(p<pi) ~ ^E^^'^ 2 ) + ^(faVM), (2-8) 

ij 

where the index i = 1,2, 3. The a 1 are the Pauli matrices, and the complex fermions if) 
obey the anti-commutation relations: 

where a = 1,2. The supersymmetry generators are now complex, but still take a form 
similar to the previous example, 

while the constraints are given by, 

C A = fABc{x B Pc ~ ^Ba^Ca)- (2.9) 

The supersymmetry algebra is now, 

{Q a ,Qp} = 0, 

{Q a ,Qp} = o, (2.10) 

{Q a ,Qp} = 25 a pH - 2a i ap x i A C A . 

The most glaring difference between this model and the nine-dimensional zero-brane case 
is that the Hilbert space is now a Fock space with a canonical vacuum. This model is quite 
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special because the Pauli matrices form a Lie algebra, and so the complex supercharge can 



be expressed as 



After introducing a potential, W = -^j 'ABC^ k x\x' B XQ, we can conjugate the supercharge 
in the following way: 

e w Q a e~ w = o^pp\. 

The study of the ground state wavefunctions then takes on a cohomological flavor since 
the supercharge acts roughly as the operator, 

Q ~ d + dWA, 

on the wavefunctions, which we can view as differential forms. We expect that, in this 
case, there should then be an explicit proof from studying the spectrum directly that shows 
there are no zero-energy L 2 wavefunctions for this model. 

The five-dimensional case is governed by the Hamiltonian: 

H = ^Tv(pY) - ^£Tr([*V] 2 ) + lTr(^ 7 V^])- (2.11) 

ij 

The index i now runs from 1 to 5, and the matrices 7 are elements of the SO (5) Clifford 
algebra. Again, the fermions are complex, and obey the relations, 

where a = 1, . . . , 4. The constraint has a form identical to the previous case (2.9), and the 
supersymmetry algebra is given by: 

{Qa,Q/3} = 0, 

{Q a ,Qp}=0i (2-12) 
{Q a ,Qp} = 28 aP H -2 1 i ap x i A C A . 
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2.3. Wrapped D-branes 



Some of the models described in the previous section have already been realized from 
wrapped D-brane configurations. Let us begin by considering type IIB string theory, and 
the case of three-branes wrapped on a collapsing three-cycle. In his study of singularities 
near conifold points of Calabi-Yau manifolds, Strominger required only a single massless 
BPS state wrapped on the vanishing cycle pi|| . That there should be no bound states 
has been argued from a somewhat different approach in [|22|j . The geometry of interest is 
Rx S 3 , where the S 3 shrinks to zero size. Clearly, the effective theory on Rx S 3 is not N=4 
Yang-Mills; such a theory would make little sense. Rather the world-volume theory of a D- 
brane on a curved space should be described by a topologically twisted theory |23|,|2| . As 
the size of the sphere shrinks, only the light degrees of freedom are relevant. The question, 
in this situation, then concerns the existence of a ground state in the theory obtained from 
the dimensional reduction of four-dimensional N=l Yang-Mills, as first mentioned in 0. 

We can also check the situation for type IIA, where we perform the same analysis for 
the case of two-branes wrapped on a vanishing two-cycle. The situation is exactly analogous 
to the case described above. The geometry is now R x S 2 , where the S 2 shrinks to zero 
size. The only difference involves the number of supersymmetries. The effective theory is 
now the reduction of N=l Yang-Mills from six dimensions. Both models were explicitly 
described in the previous subsection. It seems plausible that other D-brane configurations 
will realize many, if not all, of the remaining models which we have discussed. 

2.4- Gauge Invariance 

There is a rather nice feature of some of the computations that we shall describe that 
deserves a separate comment. Whether it provides a hint at how to formulate covariantly 
M-theory as a matrix model we leave to the judgement of the reader.! The gauge-fields 
in a quantum mechanical gauge theory are non-dynamical. They serve only to enforce 
the constraint that all states in the Hilbert space be gauge-invariant. How do we enforce 



2 When this section was originally written, the preceding comment seemed most appropriate. 
Subsequently, there has been an interesting proposal for a non-perturbative definition of the type 
IIB string, given in [ p4[ . The high temperature limit of the partition function that we describe 
in this section reduces precisely to the model in that proposal. The relation between M(atrix) 
theory and the proposal in [24] seems to be in the spirit of a 'T-duality' in the time direction. 
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such a constraint in the operator formulation? For very high temperatures, the partition 
function, 

Z(/3) = J dxtre- (3H (x,x) (2.13) 

can be well-approximated by perturbation theory. The notation that we will use through- 
out the paper may be unfamiliar, and so deserves a comment: we will often consider traces 
of some operator, say O, which we will denote as, 

trO(x,y), 

where by (x,y) we mean the usual propagation of a particle from point y to point x. In 
an explicit basis of eigenfunctions, tp n (x), with eigenvalue A n for O, this expression takes 
the familiar form, 

^2 Kipn(x)ip n (y), 

n 

where n may index a continuous parameter. 

However, in computing the partition function (2.13), it is inconvenient to try to 
trace over the gauge invariant spectrum of the Hamiltonian i.e. states \ip(x) > satisfy- 
ing Ca\iP(x) >= 0. Our first task is then to implement the projection onto gauge invariant 
states explicitly, so we can trace over the full, unconstrained spectrum. The gauge con- 
straints, Cai split into two sets of SU(N) generators: one generates rotations of the x l , 
which we shall denote C b , while the other, C^, generates rotations of the fermions. Let 
us denote the operator generating a finite gauge transformation g(t) on the fermions by 
H(g(t)) where we shall drop the explicit dependence on t. To project onto gauge invariant 
states, we insert: 



[ dt [ dxtre UACA e~ pH (x,x), 

JSU(N) J 



/ dt dxtrU(g)e~ l3H (gx,x), 

JsU(N) J 



(2.14) 



'SU(N) 

where the measure for the SU(N) integration is chosen so that j su ^ N ^ dt = 1. The trace 
is now over the full Hilbert space, including gauge-variant states. 

For small /?, we can now construct a reasonable approximation for the propagator, 

e~^(x, y) = ^L^e-^e-^e-^ + ..., (2.15) 

11 



where / = d(N 2 — 1) is the dimension of the space of scalars. We shall describe this 
approximation in somewhat more detail in the following section. The fermion projection 
operator can be expressed as U(g(t)) = e ltAC ^, which yields the expression, 

m = L W * I ^^-^e-- e -^ e .^ + .... 

As (3 — > 0, we see that the contribution from group elements away from the identity 
element is strongly suppressed. Indeed, we can then replace g by I + it ■ & , and the 
exponential term involving g becomes, 

~\it-C b x\ 2 

e 2 ? 

The term it ■ C b x is more transparent when written as ^Tr[t, x 1 } 2 , but this is precisely the 
form of a term in the potential energy, V. Indeed, in this limit, the gauge parameters com- 
bine exactly with the remaining coordinates to give a trace which is SO(d+ 1) symmetric, 
rather than SO(d) symmetric. Even the fermion projection operator combines naturally 
with Hp to give a complete symmetry between x l and t, in the computation of this trace. 
We shall put this symmetry to good use in subsequent computations. Note that for the 
case of zero-branes in type IIA, the partition function appears to arise from a manifestly 
SO (10) invariant Hamiltonian, without any hint of gauge constraints. 



3. Counting Ground States 

3.1. Defining the Index 

Ideally, to count the number of normalizable ground states for these models, we would 
like to compute the low temperature limit of the partition function, 

/dx lim tr e~ l3H (x, x). 

Except for very simple systems, that computation is beyond reach. As usual, we are then 
interested in counting the number of L 2 ground states weighted by (— 1) F where F is the 
fermion number. Therefore, we wish to compute the index, 

Ind = / dx lim tr (— l) F e~ l3H (x, x), 

J (3.1) 

= n B - n F , 
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where the trace is over the gauge invariant spectrum of the Hamiltonian. Let us first 
note that the index is perfectly well-defined. The only way that the index (3.1) could 
not be counting the net number of ground states is if the Hamiltonian had an extremely 
pathological low-energy spectrum, i.e. if the density of states diverged badly as E — > 0. 
That is certainly not the case for the models we are studying. 

Whether the index is computable is another question entirely. The purpose of this 
section is to argue that our approach to computing the index actually counts the number 
of ground states. Before discussing the issues that arise in the non-Fredholm cases, let us 
discuss in some detail the situation where there is a gap in the spectrum. First, when the 
spectrum is actually discrete, the twisted partition function is /^-independent. In these 
cases, we can compute the index in the (3^0 limit, which reduces to a perturbative 
computation. The (3 — > limit is what we will call the principal contribution to the index. 
Even in the case where the spectrum is discrete, the principal contribution, which is often 
computed as an integral over the coordinates, x, can be shifted to a boundary term. To 
see this, note that we should first perform all our analysis on a ball B R , where \x\ < R, 
and then take a limit R — > oo. We can then write, 

Ind= lim / dx lim tr (-l) F e~ 0H (x, x), 

= lim lim / dx{tv(-l) F e-P° H (x,x),+ I d(3^- tr (-l) F e -^ H (x, x)}. 
R^oo/3 o ^0j\ xl<R Jp o 6(3 

Now, computing dp of tr ( — l) F e~^ H brings down H, which we can replace by Q 2 . When 
we try to run Q around the trace, 

tr(-l) F Q 2 e-^ = -tiQ(-l) F Qe-^ H 

= -tr {-l) F Q*e-f> H - ^-tre,(-l)^e-^, 

we find that the f3 variation vanishes up to a total divergence. In this expression, e{-£-[, is 
the derivative term in the supercharge, Q. If we define e n to be the fermion in the normal 
direction to the boundary, then the index can be written as a sum of two terms, 

Ind= lim lim{/ tr (-l) F e"^ H + J / f°° d(3tr e n {-l) F Qe~^ H }. (3.2) 

i?^oo/3 o ^0 J\ X \ <R *J\x\=RJp o 

At first sight, keeping track of these various limits may seem like a technicality; however, 
that is not the case. So, rather than continue a general discussion, let us revisit an old 
friend to see how these manipulations work concretely. A number of the difficulties that 
arise in the D-particle cases will become clearer. 
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3.2. The harmonic oscillator revisited 



Let us consider a single supersymmetric harmonic oscillator, which has a unique 
ground state, and a discrete spectrum. The supercharge is given by, 

Q = ipip + ip 2 x, 

where tfjf = ip 1 = 1 and {ipi,i^>2} = 0. The Hamiltonian is one-half the square of the 
supercharge, 

= -(p 2 + x 2 - iipifo). 

Now to evaluate the principal term, we can consider the first term in (3.2), we which can 
write as, 



/. 



dx - tr^i^ 2 e~^ (x -^^ 2) + ..., 

R y/27Tf3 



where (— 1) F = iipiip2 in this case, and squares to the identity. The omitted terms are 
suppressed by powers of (3. Evaluating the trace on the fermions, or in the equivalent path- 
integral language, integrating out the fermion zero modes, gives a leading contribution in 



which gives, 



dx—=e x2 . 



If we take (3 Q to zero faster than -R -2 , this term vanishes, while if we take (3 Q to zero 
more slowly, we obtain the expected answer of one. Whether or not we get a contribution 
from this term depends on how we choose to take f3 to zero. When this term does not 
contribute, the second term in (3.2) contributes, and the principal contribution is shifted 
to a boundary term as we shall see. In this model, the principal contribution is the only 
contribution to the index. 

The boundary term gets two equal contributions from R and —R in this case, and so 
can be written, 

- / d/9tr(-^i)(^iV'2)Qe-^ 

2 J0n 
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x=R 



As R becomes large, the potential term e~^ v damps the kernel, e~^ H , for large (3. We 
therefore do not need non-perturbative information about the kernel to evaluate this con- 
tribution - a small f3 approximation suffices. Whenever there is a mass gap, we have this 
nice damping, which is the reason that the index is usually computable. 
In this case, evaluating the trace on fermions gives, 

-cso 2 



r 

\ d p x _±^ e ~^/2 



x=R 



and on rescaling, we obtain: 



/. 



OO -1 



Now we see that this term can contribute if j3 Q is taken to zero sufficiently quickly with 
R. Of course, the L 2 index is one regardless of how fast or slowly we choose to take (3 to 
zero. In the case of a model with potential V homogeneous in x with degree k, a similar 
argument can be applied. In that case, if we take j3 to zero slower than R~ k then the 
principal contribution is localized to the first term of (3.2), while if take (3 to zero faster 
than R~ k , the second term contributes. In the following section, we shall evaluate the 
principal contribution for the two D-particle case by letting (3 go to zero more slowly than 
i? -4 . This seems computationally simpler than trying to localize the contribution to the 
boundary. 

3.3. Reducing the principal term to quadrature 

To evaluate the principal contribution, we have to construct a reasonable approxima- 
tion to e~@ H . We will not need to alter the usual perturbative construction of the partition 
function because of the flat points of the potential, V. We start by writing, 

e-? H = — [ e-^-J—dz, 
2m } 1 H-z 

where 7 is a contour enclosing the spectrum of H. Let us consider the generic situation 
away from the flat points. We can approximate (H — z)~ x by a perturbation series, 

1 f e ik-(x-y) jj f 

W^- Z ^ y) = j ( k y2 + v-z) {1 ~ (k*/2 + v-z) + ■ ■ °' (3 - 3) 

where the first correction, proportional to Hp, is shown, and subsequent terms are con- 
structed iteratively in powers of (k 2 /2 + V — z)~ x . The corresponding propagator takes 
the form, 

e-^(x, y) = J^yj^e-^e-^e-^ + ..., (3.4) 
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where / = d(N 2 — 1) is the dimension of the space of scalars. This approximation is 
reasonable for small j3. The omitted terms which correct this approximation appear with 
a higher power of (k 2 /2 + V — z)~ x , in (3.3), and consequently give rise to terms suppressed 
by powers of (3 in (3.4). This approximation then suffices for evaluating the first term in 
(3.2), where we choose the take (3 to zero more slowly than i? -4 . Substituting the leading 
approximation for the propagator gives, 

lim lim / / dt tr (-l) F e _/3ff r%) (gx, x) 

R^oo f3^0 J\ X \ <R Jsu(N) 

= lim lim f [ dttv(-l) F - ^ e - ]U 2t x]2 e -f3V e -/3H Fe it A C f A + _ 

R^°°wJ\x\<rJsu(n) (2tv(3) 1 / 2 

where we have approximated the group element by I + it ■ C 6 , for reasons explained in the 
previous section. This amounts to localizing the integral over the gauge group to a small 
neighborhood of the identity - a slight twist on the usual localization to the minima of 
the potential. As usual, the inclusion of (— 1) F forces us to absorb fermion zero modes. 
In our trace, (— 1) F is realized as the volume form for the Clifford algebra (2.5). There 
are two sources for fermions: the first is from e~ f3HF , while the second source is the 
fermion projection operator, e ltAC ^, inserted into the trace. After writing the fermion 
term, UaC^ — ftHp, as ipMip for some matrix M, the trace over the Clifford factors gives 
the Pfaffian of M, which is a polynomial in x and t. On rescaling the integral, we obtain, 

j dt J dx _^ e -\^C b ^/2 e -V( X ) Ff{M ^ (3 5) 

where the integration region for t is now R-^ -1 , while for x, the region is K d ( N As 
will be clear from the subsequent explicit computation, the Pfaffian is of definite sign when 
d is odd, and the integral is thus non-vanishing. However, it is far from clear that this 
term yields an integer, and indeed, it generally is not integral. Therefore, there had better 
be a non-vanishing correction term. 

We stress again that it is very natural to consider t on equal footing with the coordi- 
nates x % . Let us denote t by x , and define 7 to be i/, where / is the identity matrix. The 
coordinates x\ now form an (A^ 2 — 1) x (d+ 1) matrix. In this notation, the matrix takes 
the form M = —(i/2)fABCX l B Y, and the integral admits an SO(d + l) symmetry which we 
shall use in section four to compute explicit values for this term in the two-particle case. 
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3.4- The non-Fredholm case 

When the Hamiltonian under consideration has continuous spectrum, the twisted 
partition function is generally /3-dependent. The heuristic reason for the /9-dependence 
is that the density of states for the bosonic and fermionic scattering states can differ. 
Supersymmetry pairs bosonic and fermionic modes, but does not necessarily preserve the 
spectral density. In these cases, the principal contribution to the index is not necessarily 
integer, and there must be an additional contribution from the second term in (3.2). In the 
case where there is a mass gap, this contribution can be perturbatively evaluated. What 
happens in the case where there is no mass gap? 

Let us choose a real supercharge, Q, which squares to the Hamiltonian up to a gauge 
transformation. For this discussion, we will set the gauge constraints to zero. Q is then 
a self-adjoint elliptic first-order operator, which anti-commutes with our Z2 involution, 
(— 1) F . Let us define Q + as the restriction of Q to the +1 eigenspace of (— 1) F , i.e. the 
bosonic states. It may be helpful to think of Q as a matrix, 



where Q* + = Q- is the restriction of Q to fermionic states. The computation that we need 
to perform is the calculation of the L 2 index of Q + . This operator, though elliptic, is not 
Fredholm. Recall that a Fredholm operator, by definition, has a finite-dimensional kernel 
and cokernel. The fact that the continuous spectrum of H — Q 2 contains scattering states 
with arbitrarily small energies implies that the image of Q+ is not closed, and the cokernel 
is infinite-dimensional, and distinct from the kernel of Q-. Hence, we take the L 2 index of 
Q + to be the dimension of Ker(<3+) H L 2 minus the dimension of Ker(Q^) fl L 2 , which is 
not, in this case, the dimension of the kernel minus the dimension of the cokernel of Q + . 

Let us consider how to compute this index. Suppose that there exists a Green's 
function, G, for Q 2 ; i.e. a self-adjoint singular integral operator G which annihilates the 
kernel of Q 2 and acts as {Q 2 )~ l on the orthogonal complement of the kernel. By definition, 
G obeys, 



where P denotes the orthogonal projection onto the kernel of Q 2 . So, P then annihilates 
all states which are not zero-energy. We recall that by singular integral operator we mean 
an operator which is obtained by integrating against a matrix- valued kernel, g(x, y), which 




Q 2 G = I - P, 
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has a well- understood singularity along the diagonal x = y. For example, the inverse of 
the Laplacian in three dimensions is the familiar kernel, ~ \x — Let us denote the 

restriction of G to the +1 eigenspace of (— 1) F by G + , where: 

G + = G(I+(-l) F )/2. 

Let P± denote the orthogonal projection onto the L 2 kernel of Q±. Then QQG + = 
(I - P){I + (-l) F )/2 = (I + (-l) F )/2 - P + , and similarly, using the fact that Q anti- 
commutes with (— 1) F and commutes with G, we have QG + Q = (I — (— l) F )/2 — P_. 
Therefore, the index of Q + can be expressed as, 

IndQ+ = trP + - trP_ 

= tr[(J + (-l) F )/2 - QQG+] - tr[(J - (-l) F )/2 - QG+Q] 

= tv((-l) F -[Q,QG+]) 

= tr((-l) F + [Q,(-l) F QG/2]). 

Of course, the difficulty is that we cannot construct G explicitly - even the claim 
that G is represented by a singular integral operator requires some justification, which we 
will give shortly. This difficulty can be summarized in the following way: given a set of 
eigenstates, ipEai with eigenvalue E under Q 2 , the inverse is formally, 

(Q 2 )~ 1= E e-^e^eM, 

Ea,E>0 

where the sum may be over continuous indices. When the continuous spectrum is not 
bounded away from zero, it is not clear that the resulting sum converges to a function in 
any reasonable sense, since G is an unbounded operator on L 2 wavefunctions in this case. 
However, if the scattering states do not pile up at low-energies, then it should be intuitively 
reasonable that G is still nice, as, for example, in the free-particle case. We will see this 
later by realizing Gas a limit of bounded singular integral operators G w . Physically, this 
limiting procedure is equivalent to adding a mass term to the propagator, and taking the 
limit where the mass vanishes. The problem we are decribing is a common one in any 
theory with massless particles. So, let us proceed along the usual path by constructing an 
explicit approximation W to G for which we can compute the trace, 

tr((-l) F + [Q,(-l) F QW72]). 
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We must then verify that for a carefully constructed W this trace is the same as the one 
computing the index of Q + . 

Our approximation will have the property that, 

Q 2 W = I-E, 

for some compact error term E given by integrating against a matrix- valued kernel, e(x, y). 
We also want e(x, y) to decay polynomially in (\x\ + \y\) to a sufficiently high power which 
we need to determine, and study in greater detail later. 

Let us describe what data is needed in order to insure that tr((— 1) F +[Q, (— l) F QW/2]) 
computes the index. This is equivalent to the vanishing of tr[Q, (— 1) F Q(G — W)/2], which 
is the difference between tr((— 1) F + [Q, (— l) F QG/2]), which computes the index, and 
tr((— 1) F + [Q, (— l) F QW/2]), which we hope computes the index. 

For any integer to, the operator W can be constructed so that the kernel for G — W 
has to continuous derivatives. Let xr be the characteristic function of a ball Br of radius 
R. The characteristic function is defined to be one on the ball and zero elsewhere. We 
will throw xr into the various traces to serve as a cut-off on the infra-red physics. By a 
similar argument to the one used in the section 3.1, we may use the divergence theorem to 
transform trx^[<5, (— 1) F Q(G— W)/2] into an integral over the boundary of Br. Therefore, 
if we can show that (— 1) F Q(G — W)(x,x) is decaying sufficiently rapidly at large \x\ = 
R, we deduce that trx#[Q, (— 1) F Q(G — W)/2] converges to zero, and the index can be 
computed by replacing G by our approximation, W. 

To prove that Q(G — W) decays sufficiently quickly, we will need to use an argument 
that may be unfamiliar to the reader which establishes a correspondence between asymp- 
totic estimates for Q 2 and decay rates for solutions ip to Q 2 ip = F, where F satisfies some 
growth constraint. In order to orient the reader, let us first examine what the argument 
says in a much simpler case. Consider the differential equation in one variable r, on (0, oo), 

where w is some constant. For this equation, a weak form of our general argument below 
says that if e awr g G L 2 , for some a G (—1, 1), then we may conclude that e awr f is normal- 
izable, if / satisfies the growth constraint e ~ bwr f G L 2 , for some b < 1. The condition on 
the growth of / is clearly necessary to rule out the addition of the non-normalizable e rw 
to any solution. In this simple case, we can prove this result by a direct integration. We 
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will use the following analogous result, which is established in much greater generality in 



Suppose that for some positive constant c and some compact set K, we have an 
estimate of the following form: 

\\Qf\\ 2 >\\cf/r\\ 2 , 

for all wavefunctions, /, which vanish outside K. With these assumptions, if Q 2 F = e, 
with r c e G L 2 , and if F satisfies the growth constraint that F/r c ~ s be normalizable for 
some positive s, then 



r 



c-1 



-F G L 



2 



l n ( r )l+e 

for all positive e. Also, 

ln(r) e 

for all positive e. In the familiar case where we have c instead of c/r, we would obtain 
exponential decay as in the one-dimensional example. The reason we have a weaker decay 
rate in this example is that the decay is roughly no worse than e - ^, where V> is a function 
with \dil)\ 2 less than the asymptotic lower bound for Q 2 , which is c 2 /r 2 in this case, and 
w 2 in the one-dimensional example. 

For purposes of illustration, let us sketch a proof of these growth estimates. For any 
function u supported in the complement of K, integration by parts yields 

0= (Q 2 F,u 2 F) = \\QuF\\ 2 -\\[Q,u]F\\ 2 . 

The assumed estimate implies that, 



\\cuF/r\\ 2 > \\[Q,u]F\\ 2 . 

From this formal inequality, we can deduce the normalizability of cuF/r when cu/r > 
\[Q, u]\. More accurately, we can deduce the integrability of (c 2 u 2 /r 2 — \ [Q, u]\ 2 )F 2 . Let us 
define a cutoff function p which is identically one outside a large ball K' containing K in 
its interior, and vanishing in K. Then taking u to be r c /ln(r) e times the cutoff function 
p gives formally, 

||cr c_1 F/ln(r) e p|| 2 < ||(r c_1 (c - e/ln(r))p + rV)/ ln(r) e F|| 2 . 
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Collecting terms, we obtain: 

J r c - 2 (2ce-e 2 /ln(r))|F| 2 /ln(r) 1+e p|| 2 <c^y |F| 2 , 

for some constant ck', which depends on K' and c. A limiting argument, approximating r 
by a sequence of bounded functions can be used to obtain from this formal inequality the 
boundedness of the left side, when |F| is bounded on compact sets. 
We will need the following variant of this inequality. Suppose that 

where w is now a positive operator with w > c 2 /r 2 . Suppose that Q 2 F = in the comple- 
ment of a compact set K. For some k to be determined, consider the point where r k F is 
maximum. By the maximum principle, we have at the maximium that, F' = —kF/r and 
(r k F)" < 0, where: 

{r k F)" =k{k - l)r k ~ 2 F + 2kr k ~ 1 F' + r k F" 

>k(k - l)r k ~ 2 F - 2k 2 r k ~ 2 F - (2d - ly^F' + r k ~ 2 c 2 F 
= (2d-k-2 + c 2 /k)kr k - 2 F. 

If we choose < k < 2(d — 1), we may deduce that if the maximum exists it must occur in 
K. If Fr a is bounded for any positive a, we may deduce that Fr h is in fact bounded by its 
values on K for all k with k < 2d + c 2 /k — 2. In our applications, c 2 will usually be the first 
or second eigenvalue of the standard spherical Laplacian, which is or (d— 1). We recall that 
the standard Laplace operator on has eigenvalues, k(d + k — 2), with k = 0, 1, 2, . . ., 
where the multiplicity of each eigenvalue is, ( 2fc +^~^K^+ d ~ 3 ) ! _ j n tn e first case, we see that 
we get r d ~ 2 decay; in the second we see that for d at least 3, we get faster than r 2 ~ 2d decay. 
These estimates extend easily to the case when Q 2 F = e, where e satisfies the condition 
that r 2d e is bounded. We will apply this to the case when Q 2 = Q 2 +Q 2 , the hamiltonian 
in the first and second variables and F will be a kernel constructed the difference between 
the Green's function and W. This argument formalizes the observation that the product 
of two elements of the kernel of H should decay twice as fast as a single element of the 
kernel, and G — W looks like such a product. 
We will show that, 

Q 2 = A r + u/r 2 , 
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acting on wavefunctions supported outside a large compact set; here r is the distance along 
the flat direction, A r is the radial part of the Euclidean Laplacian in the flat directions, 
and the operator u is semi-positive with first eigenvalue greater than the first or second 
eigenvalue of the spherical Laplacian. We can then apply the above argument to deduce 
that G — W decays like r 2 ~ 2d . We can improve this estimate by observing that we get the 
second eigenvalue when the operator is restricted to wavefunctions with odd parity, and 
the first eigenvalue for the restriction to even parity wavefunctions. Split (G — W) into its 
even and odd components. The odd component of (G — W) decays faster than r x ~ d by the 
preceding discussion, and it is easy to see that applying Q only improves the decay rate. 
To see this, use the inequality for Q 2 F = 0, 

\\wQF\\ <2\\[Q,w]F\\, 

and choose w appropriately. For the even component of G — W, we have that its image 
under Q is odd; therefore, we again have the improved decay rate determined by the 
second eigenvalue of the sphere. This estimate requires that (G — W)r a be bounded for 
some positive a and that G be given as a singular integral operator. The last condition 
is needed to ensure that (G — W) is a smooth function. We will not prove these results 
in detail but will merely sketch how they follow from the same sequence of ideas we have 
introduced. One considers first instead of G and W corresponding operators G w and W w , 
where G w is (Q 2 + -u;) -1 restricted to the orthogonal complement of the kernel of Q 2 etc. 
It is easy to get the desired initial r a boundedness for G w — Q w for each w > 0. One 
then can get the desired growth bound; i.e. we show that the max is controlled by an 
estimate on the compact set. We then allow w to tend to zero to get the desired result for 
G — W = lim w ^ (G w — W w ). This type of argument can also be used to show that G is a 
singular integral operator. 

This then establishes the desired decay estimate given one: a construction of W which 
leads to sufficiently small E; that is |e(x, y)| < + \y\)~ d ~ 1 , and two: a demonstration 
of the claimed asymptotic lower bound for Q 2 . Under these two conditions, tr((— 1) F + 
[Q, (— l) F QW/2}) computes the index. We now turn to the evaluation of [Q, (— l) F QW/2], 
which we need to boil the problem down to a concrete computation. 

It will be convenient to arrange the construction of W so that on a large compact set, 
say a ball Br, its contribution to the index can be computed by the standard principal 
term computation described in the previous subsection. On the complement of this set 
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we will need to use special coordinates to find a nice expression for W. This requires us 
to define an approximation A to QG rather than the approximation W to G, but since, 
off a compact set, A will clearly be of the form QW ', this will not affect our preceeding 
discussion. The use of two separate constructions for W in different regions may seem, 
perhaps, a bit unnatural when dealing with Euclidean space. It is forced on us, in part, 
by the need to obtain very good control of the error term E as \x\ tends to oo. This rules 
out the use of the local computation used in Br and described previously. Moreover, the 
special coordinates we use in the complement of the compact set, like polar coordinates, 
become singular at the origin. Therefore, we will need to use two sets of cutoff functions 
to patch the two approximate inverses together. 

Let p n j(x) be a sequence of cutoff functions which approach XjR,(x), and set p n = p n ,\. 
Let W be our approximate Green's function near oo, which we will construct in section 
five. The operator d(3Qe~@Q (x, y) is the standard kernel that we will use in Br. We 
can create a global approximation to QG by defining: 

rPo 2 

A(x,y) = p n (x) / dPQe-^ (x,y)p ny2 (y) + (1 - p n (x))QW'(x,y)(l - p nA/2 (y)). 
Jo 

The cutoff functions on the left of each operator are inserted to average the two operators. 
The cutoff functions on the right, however, are inserted so that the operators are localized 
to the domains where they are well-defined, and satisfy the desired estimates. The right 
cutoffs are of course chosen to be identically one on the support of the left cutoffs; otherwise, 
they would destroy the averaging effected by the left cutoffs. This is the reason for the 
second index on p. 

Then to evaluate [(-I ) F A, Q], we write: 

[(-l) F A,Q]=[Q, Pn ](-l) F f° d(3Qe-™ 2 Pn , 2 -[Q, Pn ](-l) F QW'(l- Pn>1/2 ) 

Jo 

- Pn (-1) F (I - e"^ V2 + (1 " Pn)(-l) F (I ~E)(1- p nA/2 ) 
+ (1 - p n )(-l) F Q[Q, W'}(1 - p nA/2 ) + . . . , 

=[Q, Pn ] / dpQe-™ -[Q,p n ](-1) F QW 
Jo 

- (-lf(I - p n e-^ 2 ) + (1 - p n )(-l) F (-E) 
+ (l-p n )(-l) F Q[Q,W'] + ..., 
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where the omitted terms are terms that trace to zero. We will construct W in section five 
so that the trace of (1 - p n )(—l) F (—E) + (1 - p n )(-l) F Q[Q, W'] will tend to zero as n 
tends to oo. Thus subtracting off the (— 1) F I term we are left to compute, 

tr[Q,p n ] f° d(3Qe-™ 2 - [Q,p n ](-1) F QW + (-1)V AQ '- 
Jo 

The last term is the principal term which is given by evaluating the integral (3.5). Taking 
the limit as n tends to oo, the two commutator terms converge to the boundary traces, 



/ tr e n / 

J\x\=R \ JO 



00 

d(3Qe~ PQ -txe n {-l) F QW 



Choosing j3 Q to go to zero more slowly than R 4 , the integral, 



/ tre n / 

J\x\=R Jo 



dPQe 



-PQ 
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decomposes into two pieces. One is associated with a small neighborhood of the flat re- 
gions, which consist, say, of all points of distance at most one from a flat point. The 
other contribution is associated with the complementary region, i.e. almost all of the 
sphere of radius R. It is not difficult to show that the contribution from the flat region 
is squeezed to zero as R tends to oo. The contribution from the complementary region is 
not vanishing. Using similar arguments to those presented earlier in this section, it is pos- 
sible to show without extensive computation that this term exactly cancels the principal 
term. Standard constructions give a W whose contribution from the complementary re- 
gion of tre n (— 1) F QW also exactly cancels the contribution of f^ =R tTe n d/3Qe~P® 2 . 
Therefore, the total contribution of this boundary integral to the index comes from, 

-\ I tre n (-l) F QW, (3.6) 
1 Jn f (r) 

where Np(R) is a small neighborhood of the flat points on the boundary of the space, 
which is a sphere of radius R. We will show in section five that this integral converges to 
—1/4 in the two-particle case. In summary, the additional contribution to the index from 
the boundary is computed by evaluating (3.6), which is localized to the flat directions. 
The sum of (3.6) and the usual principal term from (3.5) must then be integer. 
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4. Two-Particle Binding 

4-1. Symmetries and Pfaffians 

The simplest case to consider is N = 2; already, however, the integral (3.5) is thirty- 
dimensional for type II A zero-branes! We shall have to use the various symmetries available 
to us to simplify the computation of the principal term. Recall that the coordinates x\ 
now form a 3x (d+1) matrix. The integral (3.5) is invariant under the symmetry, x — > gxh, 
where g is an element of SO(3) acting on the left while h is an element of SO(d + 1) acting 
from the right on x. Note that the left action is a gauge transformation. 

By using these symmetries, we shall rotate x into a special form, 

61 0--- 0\ 

b 2 ••• , (4.1) 

63 • • • / 

and reduce our integral to one over only three variables. 

In these special coordinates, the potential (including the gauge parameters) takes the 
special form, 

v = - l -{b\bl + b\bl + blbl). 

We shall now evaluate Pf(M) by evaluating the determinant of M. For the moment, let 
us return to general coordinates where M = — (i/2)fABCX l B / y l - For convenience, let us 
denote x A Y by xa- The matrix then takes the form, 

• / x 3 -x 2 \ 
M = -\-x 3 Xl . 

\ X 2 -X! / 

By row manipulations, or equivalently, by studying the eigenvalue equation, we find that 
the determinant can be expressed as, 

det(M) = ^ 3( - d _ 1 ^ det(xiX2a;3) det(l — x\x 2 X x 3 x^ X x 2 x^~ x ), 

where x~ A x = (xa — 2ix A I)/\xA\ 2 - After rotating x into our convenient set of coordinates 
(4.1), we can compute this determinant, and on taking the square root obtain, 

Pf ( M ) = 22(3=1) 
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We can immediately see that when d is odd, the Pfaffian is an even function of the vari- 
ables in both special and general coordinates, and the corresponding integral (3.5) is non- 
vanishing. 

The last ingredient that we require to compute the integral (3.5) is the measure for 
our simplified coordinates (4.1). We shall obtain the measure by gauge-fixing the integral 
(3.5) using the Faddeev-Popov approach. Let us take x' to be, 




••• 

We shall insert one into the integral (3.5) in the form, 

db'dgdhS(x' - gxh)f(b), 

where we have to determine fib). For some go and ho, xo = goxho takes the form (4.1). 
The integrals over g and h then reduce to integrals in a small neighborhood of g and h , 
with the exception of the SO(d — 2) subgroup of SO(d + 1) that leaves the form (4.1) 
invariant. If T is a generator for the left SO (3) action, and R a generator of the right 
SO(d+ 1) action which does not leave (4.1) invariant, then we can replace integration over 
9, h by, 

r]{d) Yo\{SO{d - 2)) J db'dTdR5{x' -x Q - Tx Q - Rx )f(b). 
The remaining integrals are straightforward, and we find that, 

The integrals over b are constrained such that bi > b 2 > 6 3 . The symmetry factor 77(d) is 
4 for d > 2, but is 2 for d = 2 because the left and right symmetry groups are then both 
SO (3). The value of the symmetry factor can also be checked by computing a Gaussian 
integral in the <i-dimensional model, and comparing the result to the answer obtained using 
the measure for these special coordinates. 

Finally, inserting one in this form into the integral (3.5), and integrating over x,g,h 
gives for d > 2, 



vol(S0(3)) (27r) 3d / 2 



f v6L(SO(d+l))v6L(SO(3)) d _ 2 2 2 2 2 

J db 4vol(SO(d-2)) {blb2h) 1 {bl ~ h) {bl ~ h) X 
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(bi-b 2 3 )\ ¥( ^(b 1 b 2 b 3 ) d - 1 e v . 

The first factor of l/vol(5'0(3)) comes from the normalization of the integration over the 
gauge group, where we recall that we chose the normalization so that j su ^ 2 ) dt = 1 prior 
to rescaling. 

4-2. Computing the principal contribution 

Now there is a nice change of variables that will allow us to evaluate this integral. 
Set: yi = b 2 bz/V2,y 2 = bibs/ V2,y3 = bibs/y/2, and the integral becomes, 



Gfa(M 2 6 3 ) d - 2 |(^ - b 2 2 ){b'{ - bl){b' 2 2 - b'i)\ 



d-l„V 



-Ab^bsf-'e 



(27r) 3d / 2 2 2 ( d -!) 

= / dy{y iy2 y,) d -\yl - yl){yl - yl){yl - yl)^ ^ J d/2 e~ y ' +y ^ 

1 r](d-l)vo\(SO(d-3)) 



dxe -\-\ 2 l^i .... 

R3d 2 2 2 ( d -!) 7T 3d / 2 vol(50(d))vol(50(3)) 

_1 1 7/(d-l)vol(50(d-3)) 

-2 22(5=1) vol(SO(d))vol(SO(3)) ' 

Lastly, we must multiply the result by the value of tr(J) from the trace over the fermions, 
which gives an extra factor 2 3 ( d_1 ). The net result is the formula: 



P = 2 



d-2 



rj{d - 1) vol(SO(d + l))vol{SO{d - 3)) 



(4.2) 



rj(d) vol(SO(d - 2))vo\(SO(d))vo\(SO(S)) ' 

for the principal contribution, P, for d odd, where we recall that vo\(SO(n)) = 
voliS 71 ' 1 ) ■ • -vol^ 1 ), and that vol(S n ) = 27r^/r(^). Let us conclude this discussion 
by listing the explicit values in the following table: 



Dimension 


Principal Contribution 


3 


1/4 


5 


1/4 


9 


5/4 



Table I: The principal contribution to the index. 
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5. The Propagator for Well-Separated Branes 

5.1. Some general comments 

We have determined in the previous section that the principal contribution to the index 
is fractional. Since the index must be integer, there is a missing contribution. The manner 
in which this contribution arises is quite surprising, and involves a bizarre conspiracy of 
cancellations. Let us outline the procedure we will follow before presenting a detailed 
discussion. We will construct an approximation to the propagator for the two zero-branes 
when they are far apart. The approximation will be sufficiently good in the sense that 
any corrections will not contribute to the boundary term (3.6). At long distances, the 
only states that make a sizable contribution to the propagator are those localized along 
the flat directions of the potential. The simplest approximate description of the physics 
governing the light degrees of freedom is in terms of free particle propagation along the 
flat directions. This is immediately modified when we try to 'integrate' out the massive 
modes. Let us label coordinates for the d a massive directions by y. Then for example, the 
action of <9 2 , which is part of the Laplacian (2.3), on the wavefunction for the flat direction 
is modified because of its action on the harmonic oscillator ground state, 



<W 4 i ,2 /0 /r\ d */ 4 



d a | 

tt~ - \y\ 



dqjdq - 4) d a 2 \y\ A 
16r 2 4r 1 4 

Each of the terms appearing on the right is of order 1/r 2 , and so none can a priori be 
neglected. In a similar way, the rest of the terms in the Hamiltonian modify the long 
distance behavior. This includes the 0(y 4 ) terms in the potential which we recall is of the 
form, V ~ r 2 \y\ 2 /2 + 0(y ). Since we need to include the 0(y ) terms, this approximation 
is not one-loop in the usual sense. We will need to sum up all the corrections to free 
propagation, which are of order 1/r 2 , and surprisingly, they all cancel. The remaining 
index computation then involves free particle propagation on the moduli space which is 
Rjd c -2)^2 2 . From this computation, we will recover the needed corrections to the index. 

The construction that we shall describe is to be contrasted with the kind of effec- 
tive action for the light modes that has been obtained using large-distance low-velocity 
expansions [fFJTI]. After integrating out the massive modes in a one- loop approximation, 



the leading correction to the effective Lagrangian at large r is a term of order ~ i> 4 /r 7 
for the case of the nine-dimensional model. The connection between that approach, and 
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the computations that we shall describe does not seem transparent. It seems possible that 
exploring the connection in detail will give insight, and dare we hope a proof, of the desired 
non-renormalization theorem for the F 4 term. Computing the leading corrections to the 
three and five-dimensional models also seems an interesting question. Since the amount 
of supersymmetry is reduced, we might suspect that there is a correction to the metric on 
the moduli space. However, in constructing the propagator using this approach, we do not 
find any fundamental difference between the three cases. 

Let us start by discussing the form of the various operators that we need to study 
in special coordinates. Without rotating to a convenient set of coordinates, it will be 
very difficult to say anything about the structure of the partition function. We can rotate 
our coordinates, x, into a convenient basis by using a combination of gauge and flavor 
symmetries. So, we can choose a basis, 

x = kAq, 

with k G 50(3), q G SO(d) and A the following 3 x d matrix, 

It ••• 0\ 
A= y\ ••• yl . (5.1) 

\o yl ■■■ vt) 

We have set A} = r, and let y denote the remaining 2 x (d — 1) matrix, with A^ = y l A , for 
i, A > 1. The reason for this choice is that the flat directions are now at the locus, y = 0. 
Note that the choice of k and q is not unique here. So the mapping, 

m(k, A, q) — > kAq, 

projects SO (3) x R 2d_1 x SO(d) onto our space of matrices, x, but is clearly not one-to- 
one. The fibers of the map, m, are non-trivial. Any function that depends on x can then 
be lifted to a function in the product space, which constant under those transformations 
of k, A, q that leave x invariant. Most of our computations will focus on the neighborhood, 
N F , of a flat point given by \y\ 2 < 1, and r tending to oo. This is the region that contributes 
to (3.6). 

We need to write the Hamiltonian in terms of these new coordinates, which include 
the d+ 1 angular variables parametrizing the d + 2-dimensional flat directions. The kinetic 
terms in the Hamiltonian can be determined by computing the Laplacian for the metric 
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associated to this coordinate choice. Let us recall that given a metric, g, the Laplacian is 
given by: 




where \g\ is det(g), and the Tj are a basis of vector-fields for some coordinate system. First, 
we need to make a choice of basis of vector-fields. For the coordinates, (5.1), it is natural 
to have {J^, -^r}2<i<d,B>i as part of the basis. We need d + 1 additional vector- fields. 
From the left 50(3), we can choose the two vector-fields, {X 2l X 3 } which are associated 
to the two 50(3) generators, 

1 0\ / 1\ 
-10 0,0 00 
0/ \-l 0/ 

respectively. Similarly, we can add the d—1 vector-fields {V}}j>i, associated to the right 
SO(d) generators, z(j). The matrix z(j) is a d x d anti-symmetric matrix with only one 
positive entry, z(j)±j = 1. Our total basis is then composed of the subset of tangent vectors 
to the product space, 50(3) x R 2 ^ 1 x SO(d), given by {X 2 ,X 3 } U {£, ^} 2 <;<d,B>i U 
{Vjh>i. 

We need to determine the metric, g, for this coordinate choice. The set of vector-fields, 
{ Jr' df 5- } are or thonormal and orthogonal to the rest of the basis. The rest of the basis 
have inner products, 

(Xj,X k ) = r 2 5 jk + (yy^jk, 
{Vj,V k ) = r 2 5 jk + (y t y) jk , 

and, 

(Xj,V k ) = 2ry jk . 

These inner products can be determined by pushing forward the vector- fields under m, 
and computing the resulting norms. Now the metric can be written as a direct sum of two 
metrics, g = g' © g", where g" is the identity matrix for the coordinates corresponding 
to (r,y). The interesting part of the metric is the part for the angular variables. So, let 
us write, g' = r 2 I + K, where K is determined from the above inner products. Then 
(g') -1 = I/r 2 — K/r 4 + . . ., where the omitted terms are suppressed by more powers of r. 
To compute the Laplacian, we need: 
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logdet(» = logdet(#') = trlog(r 2 /) + trlog(7 + K/r 2 ) 

= 2(d + 1) log(r) + tr(K/r 2 ) - tr(K 2 /2r 4 ) + ... 
= 2(d+l)log(r)-2\y\ 2 /r 2 + ..., 

where omitted terms are again of lower order. Finally, this allows us to write down an 
expression for the Hamiltonian in these special coordinates at (k, A, q), 



d 

2H = - 



_(d+l) d_ 2yi_d__l_ v v 



j>i j>i (5.2) 



i>j>i 

where A y is the Laplacian in the y variables. We have written Hp as ^amAI^^Y^, where 
= ilapfABC^Baiicp, m the nine-dimensional model, and analogous expressions for the 
other cases. The omitted terms are all of order (9(l/r 3 ) or smaller. 

5.2. Inverting the Hamiltonian 

To invert the Hamiltonian, let us focus first on the harmonic oscillator term, H m = 
A y + r 2 \y\ 2 , with eigenvalues 2(d — 1 + n)r, where n a non-negative integer. If there were 
no cancelling fermion term, this oscillator term would immediately guarantee a potential 
linearly increasing with r, and therefore, a discrete spectrum. In order to see the cancelling 
fermion term, we write Hp = rY r + H' F , where Y r = kAiqijY^, and H' F consists simply 
of all those terms in Hp without a A\ = r factor. It is an easy task to show that the 
eigenvalues of Y r range from —2(d — 1) to 2(d — 1). In particular, when we consider 
wavefunctions for which the transverse bosons and fermions are simultaneously in their 
ground states, all terms of order r cancel, and we are left with only terms of order 1/r 2 to 
worry about in the Hamiltonian. On any state with excited oscillators, we see that H has 
lowest eigenvalue of order r. 

Our construction of W will use the nice factorization of the wavefunctions in terms 
of their behavior in the massive directions, and their behavior along the flat directions. 
We can think of H as a block 2x2 matrix, with respect to this decomposition, where the 
Hn piece corresponds to the terms in H which take the ground state for the oscillators 
back to itself. The piece H 2 2 contains terms in H which send the state with one excited 
massive boson to itself, and the off-diagonal terms act in a similar way. As a first guess, 
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one might try (and we did) to construct W as a perturbation of some nice approximation, 
Wi, to © ■ Because of the large eigenvalue of the first excited state of H m + rY r , 
almost any construction should give a good • Inverting Hn is more problematic but 
not excessively so. 

In a standard perturbative approach, we consider HWi = I — E\. Here, E\ will 
include, for example, such terms as H12W1. As a next step, set W2 = W\ + W1E1, with 
error E2 = —Ef. One could iterate this construction to construct W = W n for some large 
n, if the errors were getting significantly smaller each time. For example if each E^ had a 
kernel ek(x,x') which was bounded by (r(x) + r(x'))~ k , this would lead us to the desired 
W after some number of iterations. With this in mind, it becomes clear what terms must 
be included in our initial approximation, W±. Because any approximate will have 
upper bound of size r 2 , we see that we cannot discard any term of size bigger than or 
equal to 0(l/r 2 ) in H which either acts on the ground state, or maps an excited state to 
the ground state. For these terms lead to errors which are not decreasing under iteration. 
Actually, this is an overstatement. If a term B maps us, for example, from the lowest 
state into a higher state we see that the error will enter as H^BH^ 1 . Because H^ 1 is 
bounded above by (r(x) + r(x'))~ 1 , this term will be decreasing in r if, for example, B is 
(9(l/r 2 ). With these remarks to guide us, we must now return to analyze H, treating as 
lower order only terms which are 0(l/r 3 ) if they map the ground state into itself, or are 
0(l/r 2 ) if they mix the ground state with excited states. These lower order terms cannot 
be neglected altogether, but they simply enter as higher terms in our iterative construction 
of W. 

The gauge constraints provide some further simplication in our computations. First we 
can 'move' to the identity element k = 1 in our product space, 50(3) x R 2d_1 x SO(d) with 
coordinates (k, A, q), by a gauge rotation. From now on, we will restrict our discussion 
to the k = 1 subspace. Further, setting two of the gauge constraints to zero allows us 
to replace differentiation by Xj, which generates gauge transformations on the bosons, 
by multiplication by the fermion bilinear which generates gauge transformations on the 
fermions. For example, in the nine-dimensional case, the fermion bilinear is given by 
Qj := —T^ipisipjs- There are similar expressions for the other cases. The remaining gauge 
constraint, C\ of the three constraints in (2.2), generates a U(l) subgroup which acts on 
y. With these considerations in hand, let us turn to the task of getting rid of the massive 
modes, and obtaining an effective Hamiltonian on the flat directions which we can invert 
to compute (3.6). 
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5.3. Constructing the effective Hamiltonian 

Let us begin by computing the total contribution of terms that map the ground state 
to itself. Each of these terms will give rise to an interaction in the effective Hamiltonian of 
the form, m/r 2 , for some m, in a manner described in the beginning of this section. The 
ground state is of the form s(q)A d ~ 1 ^ 2 e~ r ^ 2 / 2 , where s(q), the fermion ground state, is 
actually a section of a bundle determined by the lowest eigenvalue of the fermion term Y r , 
since s(q) must satisfy the equation: 

{Y r (q) + 2(d-l)}s(q) = 0. 

Therefore, the ground state depends non-trivially on the right angular coordinates, q. We 
will examine the structure of the fermion ground state in some detail, shortly. Note that 
the ground state is invariant under the remaining U(l) subgroup of the gauge group. As 
a first approximation, a general state that we need to consider is a product of the ground 
state with a wavefunction, f(r,q), along the flat directions. Let us record the various 
m l /r 2 contributions to the effective Hamiltonian which acts on /, where m 3 ,m 5 and m 9 
denote the values of the contributions for the three, five and nine-dimensional models, 
respectively. First, by definition, the term H m + rY r vanishes on the ground state. 

It is convenient to rewrite y and in terms of standard annihilation and creation 
operators: 

y = 1 (a + at), 



d i t* 

^ = V2 (a ~ at) ' 

where [a,aJ] = 1. Now we can evaluate the contribution of the derivative terms in H 
acting on the bosonic ground state, 



The term in (5.2) with one radial derivative gives, 

(0|^|0) = <0| (±=± - V y ) |0> 
= 0. 
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The term with two radial derivatives then gives, 



<0|fl?|0) = 1 



d 



+ 7(°lb| 2 b| 2 |°) 



4r 2 

Let us label this contribution, m r , then rv? r = 1/4, = 1/2, and m® = 1. After noting 
that we can replace the Xf terms in (5.2) by the action of the fermion bilinears, , which 
are just matrices, we see that there are two remaining derivative terms in the Hamiltonian. 
The first is the y-^ term and it is relatively easy to analyze, 

(0|^|0) = -^(0|aat|0) 
= (l-d). 

Let us call this contribution, m y , where rriy = 1 — d. The sign of this term is critical 
since it is the only term that maps the ground state to the ground state, and gives a large 
negative contribution to the net m. 

The last derivative term comes from the action of V- 2 on the ground state, where the 
Vi generate certain flavor rotations. The (/-dependent fermion oscillator ground state is the 
only part of the full ground state wavefunction that can give a non-vanishing contribution 
in this case. So, let us now describe the ground state of the fermions in some detail. The 
real massive fermions, {V^a, V^a}, where a runs from 1 to n = 4, 8 and 16 for the three, 
five and nine-dimensional cases, can be arranged into annihilation and creation operators: 

The operators 6, obey the anti-commutation relation, b^} = 5 a p. Can we construct 
a fermion state in the kernel of Y r + 2(d — 1)? For the moment, let us pick q = 1. In this 
case, Y r = 2z7^ /3 'i/'2a' i / ; 3/3? an d we can pick 7 1 to be the diagonal element of the Clifford 
algebra: 

1 
-1 

With this choice, Y r = 2 J2^=i (^a^a — ^/ 2 + a ^n/2+a ) • Let us choose a fermionic ground 
state \0)p which satisfies, 6 a |0)F = = for a = l,...,n/2. This vacuum is 
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then in the kernel of Y r + 2(d — 1) at the point q = 1. We should point out that, when 
restricted to the (d — l)-sphere which is the SO(d) orbit of a flat point, the ground state 
takes values in a flat vector bundle, which is therefore trivial. This means that there is a 
globally defined fermion ground state wavefunction. 

To understand the fermion ground state for arbitrary q and, in particular, to under- 
stand the action of V? on the ground state, we shall study the equivalent problem of its 
action on the operator which acts as projection onto the kernel of Y r + 2{d — 1). This 
operator is given by, 

p («> = dsi>7=w (5 ' 3) 

using a contour, T, which we will take to be a small loop enclosing — 2(d — 1). This 
construction of P(q) is readily seen to be correct by diagonalizing Y r (q), and computing 
the corresponding contour integral. 

Here we should clarify that we are really interested in the action of P(q) ^ V?'-* the 
remainder is 0(l/r 2 ) and does not take the groundstate to itself. These terms will therefore 
be of interest only as perturbative corrections. The operator P(q) V? should act as a 
scalar 2m q on an appropriately chosen basis of the ground state. Our task is to compute 
the scalar. Now, P(q) is given by conjugating the projection operator at q = 1 with 
a nonconstant orthogonal matrix, whose columns give a basis of the ground state. We 
can write this as, P(q) = 0(q)P(l)0(q) t . Then V? acting on P(q) has terms where 
O(q) or 0(q) t are twice-differentiated and terms where each is differentiated once. It is 
easy to show that this last term is annihilated by P(q), and hence is not germane to this 
calculation. If the remaining term is a multiple of P(q), then the multiple will be 4m 9 , 
as both O(q) and its conjugate should contribute a factor of 2m q when differentiated. We 
preface the computation by commenting on how to compute the derivative of Y r (q). It 
is enough to compute Viqij. By a change of coordinates, we then only need to compute, 
Viqijil). This is given by: 

Vm ^ = Jt {qem)l i=o = {qV ^- 

Evaluating at q = 1 gives, 

V iqij (l) = (Vi) l3 = 5ij. 
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Now we compute: 



V — f 1 y? Y 1 dz 



where Y r is an eigenfunction of the Laplacian '^2 i V l 2 . Hence the second integrand is a 
simple double pole and therefore integrates to zero. In order to compute the remaining 
term, it is enough to make a change of coordinates equivalent to taking q = 1. Then it is 
easy to see that Vi(Y r ) takes the —2(d — 1) eigenspace of Y r to the 4 — 2(d — 1) eigenspace. 
This gives, 



£diag(^A/ i 



dz, 



(z + 2(d - 1)) (z + 2(d - 1) - 4) (z + 2(d - 1)) 

where diag(l^Y r ) 2 = P(q)(VjY r ) 2 V(q). Integrating gives - £\ diag(VjY r ) 2 /8. 

We can compute this at q = 1 where, (VjY r ) = 2iY a pip2aip3/3, and X^(^^r) 2 /8 = 
_4 7a/3' i /'2a^3/37a'/3'' i /'2a''i/'3/3'/8- Tne terms which survive P(g) are, 

-47a/?V'2a'i/'3/37a/3' i /'2a'i/'3/3/8 - 47^l/'2a^3/?7^a' i /'2/3V'3a/8 = 7a/?7a/3/ 4 

= 4(d-l). 

So, we finally get that m 3 q = 1/2, m q = 2 and m 9 q = 8. 

With the derivative terms out of the way, we can consider the two remaining operators 
in the Hamiltonian that map the ground state to the ground state. The first is the 0(y 4 ) 
term in the potential. This gives, 

<o| »ki - ylvifio) = (o|2 X>£i/£) 2 |o) 

= ^(d-i)(<*-2). 

Calling this contribution, my, we have m v = 1/4, m v = 3/2, and m v = 7. The last 
contribution comes from the kinetic term for the two gauge rotation generators, X2,X^. 
Setting the gauge constraints to zero, we can replace the angular Laplacian, Xf, by the 
operator quartic in fermions: 

n 

^2 (^l^a) 2 + (^la^3a) 2 - 
a=l 
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The only terms that map the ground state to the ground state are those proportional to 
the identity. A quick calculation gives the numerical value, n/2. Calling this contribution, 
to/, we note that mj = 1, to^ = 2 and to^ = 4. 

If we were to stop at this point, and just consider these diagonal contributions to the 
effective Hamiltonian, a quick check would show that the net to is non-vanishing. There 
would therefore be a non-trivial 1/r 2 interaction in the effective theory. Fortunately, we 
are not quite finished. First, we can shift the coefficient of the d r term in (5.2), and 
generate a new 0(l/r 2 ) term by redefining our wavefunctions. The reason this is useful is 
that on choosing an appropriate coefficient for <9 r , we can combine the radial derivatives 
with the Vi angular derivatives to obtain a Laplacian for flat space, together with some 
1/r 2 interaction. This way, we only need to deal with Euclidean coordinates, rather than 
the messier angular coordinates. We want to shift the coefficient of the d r term from d+1 
to d — 1, so we will end up with a free particle Hamiltonian on the <i-dimensional moduli 
space, together with interactions. To do so, we note that: 




where we now redefine our wavefunctions, f(r,q) = ^f(r,q). This gives us a new contri- 
bution to to, say to c = ^=-k 

So far, we have found that H acts on the ground state and first excited state in the 
following way, 

H _ (\&d + m^/r 2 b l \ 

V b H 22j' 

Here, tot = to c + to/ + my + m q + m y + m r is the total effective interaction that we have 
found so far, where m\ = 1, m\ = 4 and m\ = 16. Lower order terms in Hn have been 
omitted. Let us turn to the form of 6, which may change the effective interaction. For 
example, the terms in b acting on the ground state of order (9(r _1//2 ) cannot be neglected. 
Our initial choice of W\ is then not diagonal to the requisite order in the basis that we 
have been using. Does b contain terms of the right order? 

It is not hard to check that the only terms that map the ground state to the first excited 
state, which are of the requisite order are those involving H' F . This follows from noting that 
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H' 



' F is proportional to y which is, ~ a 



y^/r acting on the ground state. Now we see that 



we can significantly lower our energy if we define a new ground state |0)' = (I — -£p)s(q)\0). 
The factor of 2r is chosen because acting with b on a state raises its eigenvalue under 
H m + rY r by 2r. So, H\0)' = ( jA d + m^/r 2 - 6*6/2r)s(g)|0) plus lower order terms. The 
final contribution to m from b l b is computed in the same way as the other contributions, 
and we find that mf = —m^. 

Wonderfully, the 0(l/r 2 ) terms sum to zero! From these considerations, we see that 
we are reduced to a free-particle calculation. Our choice for an approximation Wu to H^ 1 
is then particularly simple: we can just take the free-particle propagator, 



5.4- Evaluating the boundary contribution 

Now that we reached the point where we have a nice simple form for Wu, we note 
that the remainder of the perturbation construction is standard. We will not belabor the 
reader with the details of this expansion, but just provide some relevant comments. As 



nice error bound. The construction of W, which is perturbative in 1/r, then follows the 
outline that we have described earlier in this section. The contribution of the Wu term is 
non-vanishing, but it is clear after extensive, arduous but standard computations, which 
involve checking powers of r, that any trace involving the rest of W will bring in more 
powers of r -1 than appear in the Wu term. These terms will therefore not contribute 
to the boundary term (3.6), in the limit where r^oo. We can now restrict ourselves to 
the free-particle Green's function, Wu- The remaining calculation is simple. We have a 
free particle propagating on R d /Z 2 . The Z 2 identification comes from the Weyl group 
action on the Cartan of the gauge group SU(2). Let us take x as coordinates for R d . 
The Z 2 action acts as parity, sending — x. It also sends the free fermions, ip la , 
where a = 1, . . . ,2{d — 1), to minus themselves. The Hilbert space of gauge-invariant 
wavefunctions is given by: 



{/o(a0, fi(x)i/> la , f2{x)ipi a ipip, ■ ■ •}• 
Each function, fa, has parity (— l) k . All we need to do is compute the boundary term, 




we observed before, any reasonable construction yields a good approximate H 2 2, with a 
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for this system. This becomes an integral over the boundary of the <i-dimensional moduli 
space: 

~\f tve n (-l) F QW 11 = uqd -L (tJ s \ tre n (-l) f ? Q-, 

2 J S d-i( r ) VOl(b a )2(d — 2) J S d-ir r \ \x — x'\ a A x'=-a 

= -1/4. 

To close this computation, let us note that the lower bound on the asymptotic be- 
havior of H is given by its lowest value on the modified lowest wavefunctions, where the 
modification involved the off-diagonal b term. The lower bound is therefore the same as 
the bound for A^, up to terms of order (9(l/r 3 ). As discussed in section three, this easily 
leads to the claimed asymptotic lower bound for H. 

To summarize: we have found a formula for the index that counts the net number 
of L 2 ground states in certain quantum mechanical systems, where the potential has flat 
directions. This involved a study of L 2 index theory for a family of non-Fredholm opera- 
tors, which allowed us to show that the prescription we presented actually computes the 
index. For the case of two-particle binding, we have shown that there is a bound state for 
coincident zero-branes in type IIA string theory. We have also found further evidence that 
there are no bound states for two-branes twice wrapped on an S 2 , and three-branes twice 
wrapped on an S 3 . Note that these models are only special points in the space of theories 
obtained by deforming the zero-brane quantum mechanics. 

The actual computation split into two parts. Computing the principal term involved 
evaluating the integral (3.5). It would be interesting, and quite non-trivial, to compute 
this integral for higher rank gauge groups. Even better would be a method for avoiding 
this integration altogether. The second part of the computation required a study of the 
propagator for the two particles when they are far apart. Surprisingly, after summing a 
variety of corrections, this computation reduced to one involving a free particle moving 
on the moduli space. Undoubtedly, there is a fundamental reason for this simplification, 
and finding it may also shed light on whether the F 4 term in the effective zero-brane 
Hamiltonian is protected from corrections. It seems likely that there will be an analogous 
reduction to a free particle calcuation for other gauge groups. As a further comment, note 
that if we had studied a system with gauge group U(l) and some charged matter, there 
would have been no boundary correction, as in the case involving H-monopoles ||. 

The sort of decay estimates that we described can probably be used to get a handle 
on the structure of the ground state wavefunction. What is needed is an upper bound on 
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how fast the wavefunction can decay along the flat directions. They may also lead to a 
vanishing theorem showing that all ground states in these systems must have a definite 
fermion number. The index would no longer be just an index, but would then count 
the total number of ground states. This would allow us to conclude that the zero-brane 
bound state is unique. Finally, systems involving marginal binding of branes with different 
dimensions can now be analyzed in much the same way. 
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